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SUMMARY 


The energy dissipated by viscosity at the edge of a 
vibrating flat plate is calculated and compared to the 
radiated acoustic energy. A correction to the Kirchhoff 
integral estimate of the noise is derived. For Helmholtz 
number of order unity and smaller the dissipation can be 
comparable to or greater than the acoustic energy. 

A viscous compressible theory of the load distribution 
on a vibrating two-dimensional body is developed. First it 
is shown that load calculations based on potential theory 
and the Neumann uniqueness condition (continuity of potential 
or pressure on the surface) are not in agreement with 
experiment or the more correct viscous theory. For a flat 
plate airfoil the eigensolution of potential theory is 
indeterminant while viscous theory yields a unique solution 
that has square root singularities at the edges. It is 
also shown that for compact surfaces the far field acoustics 
depend only on the magnitude of the eigensolutions of 
potential theory and so will be uniquely determined by the 
viscous theory. It is suggested that the general viscous 
theory of vibrating surfaces with cross-sectional geometry 
will lead to results in agreement with experimentally 
measured load distributions. 


NOMENCLATURE 


a c* (1+t) /4 


a^ speed of sound 

A see Eq. (3.58) 

b a^'U-T:)/(l+T) 


c 


chord of two-dimensional vibrating body, 
see Fig. 2 


Co(t),Ci(t) time dependent coefficients of eigensolutions 
of flat plate problem, see Eq. (3.28) 

AdB see Eq. (2.28) 

energy radiated acoustically, see Eq. (2.23) 

E viscous dissipation energy, see Eqs. (2.23) 

and (2.21) 


F 

F(t) 


G(r) 

h 


4 ^^ 

-h 

1, J ,k 


I 

k 


force on a vibrating flat plate, see Eq. (3.29) 

viscous correlation function, see Eq. (2.10) 
and Appendix A 

Greens functions for the viscous compressible 
problem, see Eqs. (3.64) and (3.65) 

acoustic correlation function, see Eq. (2.24) 
and Appendix B 

thickness of a vibrating elliptic body, 
see Fig. 2 

perturbation enthalpy 

Hankel function, see Ref. 4 

unit vectors along the x,y,z coordinate 
axes, see Fig. 1 

see Eq. (2.14) 

co/ a 


modified Bessel function, see Ref. 4 
a ' (x) complex surface loading 
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^dB 

see Eq. (3.58) 


M 

moment on a vibrating flat plate, see Eq. (3. 

n 

unit normal to S (or C with 

zero subscript) 

p' 

perturbation pressure, (= p^h') 

q 

time average of quantity q 


Iql 

absolute value of quantity 

q 

Q 

density of viscous dissipation energy, 
see Eq. (2.2) 

r.e 

polar coordinates 


s 

see Eqs . (3.36) and (3.37) 


Sf 

far field spherical surface 
the flux of acoustic energy 

for calculating 

t 

time 


? 

unit vector tangent to C , 

see Fig. 2 

Tn(X) 

Chebyshev polynomial of the 
see Ref. 4 

first kind. 

U^(X) 

Chebyshev polynomial of the 
see Ref. 4 

second kind, 

U(s) ,W(s) 

surface velocity components of a vibrating 
body, see Eqs. (3.61) and (3.62) 

-A 

perturbation normal velocity on 

V 

domain bounded by the vibrating surface 
and 

W 

1 

CD 

0) 

CL 


Wo(t) 

vertical translation velocity 

w^(t) 

amplitude of pitch vibration 

W 

rate at which work is done 
by a vibrating surface, Eq. 

on a fluid medium 
(2.1) 

x,y 

vector field points 
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X,Y 

z 

a 

T 

<S(x) 

K 

X 

V 


Pco 

T 

<P 

$ 

0) 

S' 

div 

grad 

Im 

Re 

(*) 


I <^«y) 

X + iy 

see Eq. (3.46) 

Euler's constant, = 0.57721 
Dirac delta function 

effective Helmholtz munber, see Eq. (2.31) 
/o)/2v 

kinematic viscosity 
density of fluid medium 
h/ c 

velocity potential 
complex potential 
see Eq. (2.4) 
vibration frequency 
perturbation vorticity 
vector divergence operator 
vector gradient operator 

denotes imaginary part of a complex quantity 
denotes real part of a complex quantity 
Laplace operator 
conjugate of a complex quantity 
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I . INTRODUCTION 


The results reported in this dociament supplement those 
reported in Refs. 1 and 2. The overall objective is to 
incorporate viscous effects in the calculation of noise 
produced by a thin surface when it vibrates or interacts with 
an xonsteady flow or an incident acoustic field. The results 
are of two types. First, we derive an explicit formula for 
the viscous dissipation at the edge of a vibrating three- 
dimensional surface with constant load and give an estimate 
for the excess noise that would be predicted with the 
Kirchhoff integral. The results are a direct extension of 
the two-dimensional results reported in Ref. 2. 

Second, we investigate the surface load problem in 
detail and its relevance to the acoustic problem. The 
importance of the Neumann boundary condition and its 
similarity to the Kutta condition as a uniqueness criterion 
is illuminated. The invalidation of the Neumann condition 
is demonstrated with the surface loads data measured by 
Brooks (Ref. 3). The correct viscous origin of the load 
is demonstrated for a vibrating flat plate. Finally, a 
detailed formulation of the viscous compressible 
surface loads problem for a surface of finite thickness is 
presented. 

II. VISCOUS DISSIPATION AND THE KIRCHHOFF CORRECTION 
FOR THREE-DIMENSIONAL SURFACES 

The following derivation is the three-dimensional 
counterpart of the two-dimensional results derived previously 
in Ref. 2. The end result is a simple formula for correcting 
the noise estimate that would be obtained with a Kirchhoff 
type integral . 

Consider a vibrating three-dimensional surface S in a, 
stationary acoustic medium (see Fig. 1). From Ref. 1, Eq. 
(2.13), the time averaged rate at which work is done on the 
fluid medixim by the vibrating surface is 

W = J' h'v^ dz + Q dV (2.1) 


where h' is the perturbation enthalpy and v^ is the 
perturbation velocity on the far- field surface S^ . The 
first term is the energy radiated to the far field in the 
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Point P 



Figure 1 - Three-Dimensional Vibrating Surface 

form of sound. The volume integral is the energy converted 
irreversibly into vorticity by viscosity. We neglect 
compressible viscous effects so that the time averaged 
dissipation can be expressed in terms of the square of the 
vorticity, i. e. , 


Q = v|w' l'^' (2.2) 

If 5- ' (y) is the complex surface loading (force/unit area) 
then w' can be expressed as follows: 


to' = S X grad ReCij^e^'^^) 


(2.3) 


where 
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ip = 


1 



(2.4) 


and 0 ) is the vibration frequency, is the density of the 

fluid mediiom and S denotes the vibrating surface (see 
Fig. 1). 

Substitute Eqs . (2.2), (2.3) and (2.4) into the second 
term of Eq. (2,1) to obtain the dissipation energy; i.e.. 




/ 


dV 


5 — J y i'(y)dy y «'*(y')dy' 


( 2 . 6 ) 


y^ 


x-y 


a I x-y 


V 


x-y 


dx 


(2.7) 


I x-y 


where 


,2 . + 3 ^ 


9y^ 972 


y = (71.72) 


( 2 . 8 ) 


and the asterisk denotes complex conjugate. The volume 
integral is evaluated in Appendix A so that 


7 



a 


E 


V 


Sirp^w 


y* ^''(y)dy J S.'*(y')dy' 

S S 


•V^*F(X |y-y* |) 


(2.9) 


where 


and 


F(t) = ^ 


sin 

t 


t 


t > 0 


( 2 . 10 ) 



( 2 . 11 ) 


The viscous dissipation has thus been reduced to a 
double integral over the surface loading. The co rrela tion 
between load points on the viscous length scale /2v/o) is 
the origin of the dissipation. Also, we will show that the 
principal contribution to the dissipation arises from the 
edge loading. 

Consider the special case of a constant load surface; 

i. e. , 


I' ~ ~ Constant 


( 2 . 12 ) 


Then 


E 


V 



where 


(2.13) 


I 


y dy' y* dy V^F(X|y-y'|) 
S S 


(2.14) 


The second integral can be reduced to a line integral around 
the contour C ; i . e . , 
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I 


(2.15) 


= y* dy' • ds F' (X ly^-y' |) 
S C 


Myp-y' l-Hp 

lyp-y' 1 


where the zero subscript denotes a quantity evaluated on 
C . Now interchange the surface and contour integrals and 

let y' = yp "h z . The result is 


I 





dz 


(2,16) 


For a >> 1 , the support of the viscous correlation function 
F is of order 1/a so that the surface integral can be 
approximated as follows. Assime that the contour C is 
smooth and introduce polar coordinates at a typical point on 
the contour as shown in the inset in Fig. 1; i.e., 


z = (r COS0, r sin0) 


(2.17) 


Then 


3tt/2, 


I = 


y* ds y* rdr J' de F'(xr)(-x cos0) 

tt/ 2 


C o 


= 2 


ds / 
C o 


Xr F' (Xr) dr 


- I / F'(t)dt 


(2.18) 


C o 
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But 


/ 


t F' (t)dt = - f 
o 


F(t) dt 


so that 


o 


I 




(2.19) 


( 2 . 20 ) 


where C is the length of the contour C . Finally substitute 
Eq. (2.20) into Eq. (2.13) to obtain the dissipation, 


E 


V 


2 

o 


16 




Constant Load 
Surface 


( 2 . 21 ) 


Note that the three-dimensional dissipation is independent 
of the viscosity coefficient as we found previously in the two- 
dimensional case. (See Eq. (40) of Ref. 2.) 

While Eq . (2.21) was derived rigorously for a constant load 
surface it is not difficult to generalize the result for a 
load dis trib ution that varies slowly on the viscous length 
scale /2v/o) . The more general result is 


E 


V 


1 

16 p^o) 



Slowly Varying 
Load 


( 2 . 22 ) 


This result will now be used to derive a correction to the 
Kirchhoff integral estimate of the noise radiated by a 
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vibrating surface. From Appendix B, Eq. (B.30), the acoustic 
energy radiated by a vibrating surface is 


f 

' (y) dy 

^a ., 22 ^ J ^ 

32it p„a^ g 

where 

2tt 

r 

tt/2 

G(r) =1 d(J) 

kJ 

J 

o 

-T\ll 


j Jl'*(y')dy'- 


G(y-y’) (2.23) 


2 ik COS0 X *r (2,24) 

desin^6cos6e 


and 


k = w/a. 


x„ = 1COS0 + 1 sin0 
s 


(2.25) 


Note the similarity of Eqs. (2.9) and (2.23). The acoustic 
energy is a double integral over the surface load with a 
correlation function G on the length scale ao/w . The 
viscous dissipation is a double integral over the surface 

2 

l oad w ith a correlation function V F on the length scale 
/2v/w . For compact surfaces the two integrals are comparable 
in magnitude. To illustrate this important point we consider 
an acoustically compact surface in which case (see Appendix B) 


E 


a 


24 TTp^a 



2 


(2.26) 


Consider the quantities 



(2.27) 
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and 


AdB — - 10 log^Q ^ 


(2.28) 


which is the dB increment by which the noise will be 
estimated by a Kirchhoff formula. Then 



E 



E^/Ev 


and 


k 


d^! 


24TTp^-ao 


(16p“k-a^) 




= K*- 


£ 1 ^ds 

o ' 


where 


Finally 


and 



AdB = - 10 login 

1 + K 


over- 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


(2.33) 
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The Kirchhoff noise overestimate depends on the third power 
of the Helmholtz number based on a length scale defined by 
the square bracketed quantity in Eq. (2.31). For compact 
surfaces this scale only depends on the surface and edge 
loading and the geometry of the surface. It is clear from 
Eq. (2.31) that a fixed area surface with a very large 
perimeter will generate a smaller effective Helmholtz number 
and thus more dissipation than one with a minimxam perimeter; 
e.g., a circular disc. Also it is evident that the Kirchhoff 
noise correction factor becomes small as the effective 
Helmholtz ntjmber exceeds unity. The correction is most 
important in the compact regime. The formula (2.31) can be 
used to correct Kirchhoff noise estimates for arbitrary 
compact surfaces and generalizations to moderately non- compact 
surfaces can be derived with the more general expression, 

Eq. (2.23), for the acoustic energy. 

It is not our objective here to derive a collection of 
specific formulae with the foregoing results. Rather, it is 
our intention to shed some light on the deeper question of 
how to calculate the surface and edge loading. In reality 
the edge loading must go to zero on some scale. However, 
we have seen from the Brooks experiment. Ref. 3, that this 
scale must be very small. The edge load to be used in the 
preceding formula can be estimated from experiments like 
those of Brooks. However, this would be a costly endeavor 
if we must resort to experiment for every case and it 
becomes highly desirable to derive a theory that will enable 
us to calculate the edge load. An important conclusion based 
on the viscous and acoustic energy estimates in this report 
and Ref. 3, and the results of the Brooks experiment (Ref. 3) 
is that a correct load distribution theory must account for 
viscosity. That is the subject of the following section. 

III. THE SURFACE LOAD PROBLEM 

A. A Two-Dimensional Potential Problem 
and the Neumann Condition 

Consider the typical elliptic two-dimensional body that 
vibrates along the vertical or y-axis as shown in Fig. 2. 

We consider the fluid medium to be inviscid- incompressible 
and calculate the surface loading from the solution of the 
following potential problem: 



(3.1) 
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Figure 2 - Vibrating Elliptic Cross Section 


34 ) 

9n 



Wo(t)j*n 


on C 


4> continuous on C (Neumann Condition) 
with the perturbation pressure given by 


P 


( 


34 ) 

3t 


(3.2) 


(3.3) 


(3.4) 


The purpose of this section is to investigate the consequences 
of the Neumann Condition (3.3) as the elliptic cross section 
tends to a flat plate. 

To solve the boundary value problem we consider the 
complex potential 


$(z)=4) + i4) , z = x + iy 


(3.5) 


and let 


z 


w + b'^/w 


w 



(3.6) 


map the body into a circle 


w 


ae 


.10 


(3.7) 
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in the w-plane. The surface of the body is given by 


i0 . b^ -i0 
z = ae + — e 
a 




a + — jcosO + i (a - — ) sin6 
a I \ a 


\ . 

I SX1 


(3.8) 


so that 


a + b^/a = c/2 


a - b^/a = h/2 


(3.9) 


or 


a = c 


1 + T 

~^r~ 


b = a 


JTZI 

V 1 + T 


(3.10) 


with 


T = h/c 


(3.11) 


In the w-plane, it is easily shown that the problem 
for $ is 


= 0 


(3.12) 


Re 


3p I 

k / n = 


1 + T ”0 


w (t)sin 0 


(3.13) 


p=a 


# continuous on C 


(3.14) 
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and the solution is 


2 i 

^ = - - (1+T)W 


(3.15) 


With Eq. (3.4) the pressure on the surface is 


P' 


2 Pco^ 

(1+t) 


w^(t)sin0 


(3.16) 


But from Eq. (3.8) 


X 



COS0 


(3.17) 


and 


sin 0 




for 0 < 0 < 7T 


(3.18) 


Thus on the upper surface 


P 


2p a 
1+T 


W^(t) 




2 


or using Eq. (3,10) 


P' 



Wo(t) 



(3.19) 


The main points to be made about the solution are: 

1) It is unique. 

2) It is continuous on C. 

3) Even for the flat plate the pressure tends to 
zero elliptically at the edges. 
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I 


All of the above results are due entirely to the requirement 
of the Neumann condition; i.e., that the potential (or 
pressure in the present case) be continuous on the boundary. 

In the next section we consider an alternative approach 
to the flat plate problem in which the Neumann condition is 
relaxed. We should point out here that the Neumann condition 
is strictly a mathematical requirement much analogous to the 
Kutta condition that insures uniqueness of the boundary 
value problem. The physical basis or in fact the invalidation 
of this condition will be considered in Section III-C. 

B. Vibrating Flat Plate - Inviscid 
Incompressible Fluid Medium 

Consider the flat plate illustrated in Fig. 3. We suppose 
that the plate can vibrate in vertical translation with velocity 
WQ(t) or pitch about mid chord with angular velocity Wj^(t). 



Figure 3 - Vibrating Flat Plate 
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The boundary value problem is 



( 3 . 20 ) 


9(j) 

ay 


y=0± 


W^Ct) + Wj^(t) — 


(3.21) 


grad 4> ~ 0 at °° 


(3.22) 


No requirement of continuity of cj) on C 

The solution of Eq. (3.20) that satisfies Eq. (3.22) can be 
expressed as a superposition of dipoles on C ; i.e., 


where 


c/2 

* - - ^ J A4. in R d5 
-c/2 


(3.23) 


Alp = ip(x,0~) - (p(x,0'^) 

R = l(x-C)^ + (3.24) 

Now apply the boundary condition Eq. (3.21) to obtain the 
integral equation 

1 

wk J ‘ - ‘^<”0 + 

-1 
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where the "slash" indicates a Cauchy principal value integral 
and X = 2x/c is a normalized coordinate. Integrate Eq. (3.25) 
with respect to X to obtain 

1 

^ dY - - I ^ U2<X) 

-1 ^ 

where Chebyshev polynomial of the second kind 

and Cj^Ct) is an arbitrary constant. The general solution of 
Eq. (3.26) is well known (e.g., see Ref. 1, p. 28); i.e., 



+ C^-U^ (3.26) 


A(j) 


Co + - c/2(WqT2(X) + W3^/4 T3(X)) 


(3.27) 


where C^(t) is a second arbitrary constant that arises 
from the eigensolution of Eq. (3.26), 

The pressure on the upper surface of C is given by 

I 9Acj) 

P - - Pco Tt T ~JT 


Poo 

T 


+ c 

O 


Cj^Tj^(X) - c/2(WqT2(X) + W3^/4 T3(X)) 


y 


l-X^ 




(3.28) 


For arbitrary values of the coefficients *^o’^l pressure 

has a square root singularity at both edges of the plate. 

Also we note that the total force and moment on the plate are 

proportional to and C^ respectively. That is, for the 

vertical force 
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c/2 1 

F = I Ap ' dx = - c I p ' dX 

-c/2 -1 


I 


— " c C 


I 


dX 


\/l-X^ 


00 ✓ X 

= - -j- • TTC • C^Ct) 


(3.29) 


and for the moment about x = 0 , 


c/2 

M = f xAp ' dx = 


/ 

■c/2 




Xp' dx 


-1 




- - ~ 2 ~ ’ (" 2 ^ * C-|^(t) 


(3.30) 


The existence of an eigensolution of the flat plate 
problem is a direct conseqence of the relaxation of the 
Neiamann boundary condition; i.e., continuity of on the 

plate. To illustrate this point we apply the Nexamann condition 
to Eq. (3.28). for continuity of p' at x = ± 1, we must have 
the numerator of Eq. (3.28) vanish at both edges; i.e., 



+ C3_T^(±1) 


c 

I 










0 
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or 


• “2P (3.31) 

For pure translation we have w^ = 0 and so 


• vA (Neumann \ ,o on\ 

P -^""oVl-X ^Conditionj 


which result is the same as Eq, (3,19), In the previous 
section we built the Neumann condition into the solution by- 
choosing the circle plane solution as the standard. There 
is no a priori reason to choose the Neumann solution. In 
fact we will show in the next section that the solution of 
the viscous problem always has singularities at the edges. 
For mathematically sharp edges, the singularities weaken 
asymptotically as the viscosity tends to zero. 

Before we attack the viscous problem, we show one 
point about the compact acoustic problem that uses the 
last result. The acoustic pressure is given by 


where 



R = 


[(x-?)^ + 


o 1/2 

y ] 


(3.33) 


(3.34) 


k = to/a^ (3 . 35) 

and (z) is the Hankel function (see Ref. 4, Chapter 9). 

In the far field we get 
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Pa ■ ^ ^Poo Vife e 


y; 


■i (kr-TT/4) 


(3.36) 


where 


c/2 

S = J A* 


c/2 


1 

I 

-1 


^ * A(f) dX + ik(^) cos 9 J' A(f)*X dX 

-1 


or 


S = 


J • ttCq + (j) Ik 2 Cj^cose 


(3.37) 


and 


» r:: 


TTckp^ — -i(kr-TT/4) 


Pa 5 » Trier 


sin 6 


(i(jjC^) + i ( icjjC j^) CO s ^ 


(3.38) 


for a compact plate. Thus the acoustic far field is completely 
determined by the magnitude of the eigensolution and there- 
fore ultimately by the "uniqueness criterion" that we impose 
to solve the surface load problem. The Neiimann condition 
is one alternative but in fact is only a mathematical condition 
that has no physical basis. The load distribution measured 
by Brooks (Ref. 3) is a prime example of the invalidation of 
the Neumann condition. Any viable alternative to the Neumann 
condition must account for viscosity in some rational way. 

That is the subject of the following section. 
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C. Vibrating Flat Plate - Viscous Compressible 

Fluid Medium 

To further explore the importance of the eigensolution 
of the inviscid problem we turn now to the more general problem 
of a two-dimensional flat plate that vibrates in a viscous 
compressible fluid medium. The boundary value problem is 
stated as follows. 


1 9h ' _i_ j • "^1 A 

^ + div V =0 

a 

o 


(3.39) 


8v' 


^ + grad h' = - V grad ' x k 


(3.40) 


J -v 


I 

!• V 


= w(x) e 
= 0 


iwt 


on y = O' 


(3.41) 


Outgoing or damped waves at “ 


(3.42) 


All dependent variables are proportional to 


e^^^ so that 


V^h' + k^h' = 0 
- a^Q' = 0 


where 


k = o)/a 

o 


a 



iiT/4 

e 


(3.43) 

(3.44) 

(3.45) 

(3.46) 
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and is the perturbation vorticity. It is easily shown 

that the solutions of (3.43) and (3.44) that satisfy the 
field equations, the no-slip boundary condition and the far 
field boundary condition are of the following forms : 


c/2 

h' = - I J Ah' ^ H^2) 
-c/2 


(3.47) 


Q' 



c/2 

j 4 K„(aR)« 

-c/2 


(3.48) 


Now apply the boundary condition on the normal velocity 
component at the surface to obtain the following integral 
equation : 


c/2 

^ / 

-c/2 


AhV 


8x L 


H^^^(k|x-C|) + KQ(a|x-?|) 


T % 


-1 


^ kV2)(k|x-5|) 


dS 


= - 2iw W(x) 


(3.49) 


We demonstrate next that Eq. (3.49) has a unique solution that 
is singular at both edges. 

Consider the plate to be acoustically compact and of 
sufficiently small chord that 
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« 1 


(3.50) 



Then the kernel of Eq. (3.49) can be reduced asymptotically to 
obtain the following integral equation in terms of normalized 
coordinates? 


1 

TT 


1 

I 

-1 


Ah' ilnlX-YI + £n 


ac 

"zr 


+ Y 


dY = 


8v 


W(X) 


(3.51) 


where y is Euler's constant (= 0.57721) and 


(x,y) = f (X.Y) 


(3.52) 


Differentiate Eq . (3.51) with respect to X to obtain the 
familiar airfoil equation 


1 

dY = -^ W (X) (3.53) 

-1 

whose general solution is 



Ah 


, - ^ 


8v/c 


n/i-x2 -1 


A 


-Y^ W’ (Y)dY 


X-Y 


(3.54) 


whe re is an arbitrary constant. Now multiply Eq. (3.51) by 

1/ ^ 1-X^ and integrate over the plate. The result is 
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1 AT-, * riY = - 

8v/ TTC 


i 

r W dY 


1 lALl Cll — 

-1 

( to ^ + Y 

-1) 

J /l-Y^ 

(3.55) 


1 


which serves as a natural 
A„ : i . e . , 


integral constraint to evaluate 


1 

IT 


/ 

-1 


Ah'dY = A 

o 


(3.56) 


The solution of the viscous compressible flat plate problem 
is 


Ah' = 

\/ 1-X^ 



I 

-1 


1-Y^ W (Y) 
X-Y 


'N 

dY ► 




(3.57) 


Remarks': 

1) The solution of the viscous compressible vibrating flat 
plate problem is unique. 

2) The load distribution is of the same form as the eigen- 
solution of the inviscid problem. Furthermore the load 
has square root singularities at both edges. 

3) The vorticity is singular at the plate edges. 

While the foregoing remarks are rigorously proved in the 
asymptotic regime kc/2 << 1 and |a|c/2 << 1 the results remain 
true in the general case. The reason is that the kernel in 
Eq. (3.49) has a logarithmic singularity when viscosity 
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is included. The general theory of such equations (Ref. 5, 
for example) proves the uniqueness of solution and the 
existence of edge singularities. 

In Fig. 4 we illustrate the numerical solution of 
Eq. (3.49) for k = 0 and |a|c/2 = 4. Even for this 
relatively large value of the viscous parameter the load 
distributions for a plate vibrating in translation or pitch 
have singularities at the edges. The results in Fig. 5 are 
for a plate vibrating in pitch with |a|c/2 = 10. The 
singularity in the real part of the load is starting to weaken 
while the imaginary part still remains singular. Finally in 
Fig. 6 we show the results for translation when |a|c/2 = 100. 

The real part of the load has the vl-X^ dependence that is 
obtained with potential theory and the Neumann condition (see 
Section III -A) . The imaginary part of the load tends to zero 
as the viscosity becomes smaller but appears to have a singu- 
larity at the edges. 

The main conclusion of the viscous theory is that the 
physical effect of viscosity is absolutely essential to 
formulate the surface loads problem correctly; i.e., so that 
a unique solution can be obtained without a mathematical 
artiface. The second conclusion is that for the case of 
mathematically sharp edges, viscous effects will not eliminate 
the edge singularities. The effect of edge geometry and/or 
non-linear fluid effects must be taken into account to obtain 
a singularity free solution. Another implication of the 
viscous theory is that the Neiamann condition is not the correct 
physical uniqueness criterion to apply to the potential theory 
of thin vibrating surfaces. A strong piece of experimental 
evidence that also supports this conclusion was obtained by 
Brooks (Ref. 3). The measured surface load distribution for 
an airfoil vibrating in pitch about the 407o chord location 
is illustrated in Fig. 7. Note the relatively flat nature 
of the loading near the edges. In fact, the load appears to 
increase slightly near the leading edge. The experimental 
results indicate the existence of a strong vorticity field 
near the edges and negates the validity of the Netimann 
condition. Recall Ref. 2 that we used this evidence to explain 
the discrepancy between the measured and calculated (with the 
Kirchhoff integral) far field noise. 

We conclude this section with a simple result based on 
the potential theory of Section III-B. Consider the normalized 
load distribution for a plate vibrating in pure translation 
(Eq. (3.28)) and expressed in dB ; i.e.. 
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FiguiTG 4 - Complex Suirface Lioad. on a Flat Plate Vibnating 
in Translation and Pitch - Viscous Theory 
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Re 



Figure 6 - Complex Surface Load on a Flat Plate Vibrating 
in Translation - Viscous Theory 
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5PL dB(ref. 20 u Pa) PHASE <l>p - <t>a. deg 



Figure 7 - Surface Load on a Vibrating NACA 0012 Airfoil 
(Measured by Brooks, Ref. 3) 
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Edge Load dB Arbitrary Scale 


^dB “ 


(3.58) 



where is an arbitrary reference pressure. When A = 0 , 

the Neumann load distribution is obtained. We compare the 
Nexamann result with the case A = 0.4 in Fig. 8. The 40% 
increase in the magnitude of the Neumann eigensolution yields a 
load distribution that is very flat near the edges. By con- 
struction it does indeed have square root edge singularities. 

The important point is that by giving up the Neumann condition 
as a uniqueness criterion we can fit experimental surface loads 
data much more accurately. We should point out that this point 
is also true in regard to the Kutta condition when applied to 
the potential theory of lift. By giving up the Kutta condition 
Pinkerton (reported in Ref. 6) found that much better agreement 
could be obtained between measured and theoretical surface 
loads data. More recently, it has been shown (Ref. 7) that 
by including viscosity and section geometry (real and effective) 
in thin airfoil theory that both steady and unsteady 


A = 0.4 



-.8 X --6 


Figure 8 - Comparison of Edge Loading with and Without 
the Neumann Condition (See Eq , (3.58)) 
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results can be calculated uniquely and accurately. It is the 
fundamental contention of the present study that viscous and 
geometric effects must be included simultaneously in a 
complete theory of the surface load required for acoustic 
calculations . 

D. Vibrating Two-Dimensional Body - Viscous 
Compressible Fluid Medium 

We conclude our present study with a precise formulation 
of the viscous compressible surface load problem for an 
arbitrary two-dimensional section. The boundary value problem 
is very similar to one posed in Section III-C, Eqs. (3.39) 
through (3.40). Referring to Fig. 9 we have 


1 


ah' 

at 


+ div v' = 0 


(3.59) 


av' 

at 


+ grad h' = - V grad f2'x 




(3.60) 


ioat 

n*v'=W(s)e on 


(3.61) 
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(3.62) 


lo)t 

t*v'=U(s)e on 


Outgoing or damped waves at infinity (3.63) 


Equations (3.43) through (3.46) are still valid. We introduce 
two Greens functions G„ and G„ as follows ; 


V^Gjj + ~ 2776(x-y) 


(3.64) 


V^G^ - = 2iT6(x-y) (3.65) 

The solutions are required to satisfy the boundary condition 
Eq. (3.63). Now use the following identities: 

Gj^V^h' - h’V^Gj^ = V(Gj^Vh' - h'VGjj) 

= - 27rh'6(x-y) (3.66) 

and 

G^V^J2' - ^2'V^G^ = V(G^7Q' - 

= - 2nQ’S(x-y) (3.67) 


Integrate the above identities over the domain exterior to the 
body to obtain the following representation of the solution: 
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Next we project the momentum equation onto the surface 
to obtain the following pair of equations: 

l„U(s) + - V on (3.70) 

lioW(s) + ^ - - V ^ on (3.71) 


(3.68) 


(3.69) 


These equations can be used to eliminate the normal derivatives 
of h' and ' in Eqs. (3.68) and (3.69). After integrating 
by parts in each term that involves a tangential derivative, we 
obtain 


h'(x) 


= A. ^ 

2^ J 


ds 


9 Gy 9 Gy 

^ ^ - i‘^GHW(s') 


I 


(3.72) 


vf2' (x) = 




9G„ 

j^+va' l«.GgU(s') 


] 


(3.73) 


The final step is to take the limit as x ->■ Xp(s) - any point 
on the body contour . Note the following results: 
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2it 


ds' 


2it 

y -^0 


1 

“ 1 


(3.74) 


Lim ^ 
^ 2tt 

X^X 


s4*£ 

/ 


s-e 


3G 

air 


/ 

-a 


tW 


dt = 


where G can be or G^ since both functions behave 

like £nlx| near the origin. With Eq. (3.74) it is easily 
shown that 


ds’Gj^W(s') 

(3.75) 

and 

va'(s) - - ds’ (h’ VI!' M^ds'Gj,U(s') 

(3.76) 



This is a coupled pair of inhomogeneous singular integral 
equations for the unknown pressure and vorticity. The combined 
effects of viscosity and geometric cross section are 
rigorously included in the formulation. It is clear from our 
previous discussion in Section III-C that the solution of 
this pair is unique. If the contour is smooth, no singularities 
will appear in the solution. It is firmly believed that the 
solution of this pair of equations has the strongest possi- 
bility, of yielding results in agreement with the load 
distributions measured by Brooks (Ref. 3). The only physical 
effect that we have omitted is the possibility of nonlinear 
flow near the edges. Since there is nothing in the experi- 
mental results of Brooks that suggests the slightest bit of non- 
linearity, we feel that such effects can be safely omitted. It is 
recommended that the numerical solution of the pair of 
integral equations derived herein be the subject of a future 
research effort. 


IV. CONCLUSIONS 

The main conclusions of this study are sijmmarized below: 
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1. For Helmholtz numbers of order unity or less the energy 
dissipated by viscosity at the edges of a vibrating 
three-dimensional plate can be comparable to or greater 
than the acoustic energy radiated to the far field. 

2. The Neumann boundary condition is not the correct 
uniqueness criterion to apply to the potential theory 
of the surface load distribution on a thin vibrating 
plate. The viscous theory and the experimental results 
of Brooks (Ref. 3) support this conclusion. 

3. The correct formulation (to obtain a unique solution) 
of the surface loads problem must include the physical 
effect of viscosity and the effect of geometric shape. 

4. The far field acoustic solution is determined uniquely 
by solution of the viscous surface loads problem. 
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APPENDIX A 


VISCOUS CORRELATION FUNCTION 


Consider the integral 


T(x-y)T (x-y )dx 


where 


T(x) = ^ 


-A 1x1 


x2 = i“ 

V 


We introduce the three-dimensional Fourier Transform pair; 
i. e. , 


r 

I -la* 

.a) =J e 


^ T(x)dx 


T(x) = 


f ia- 

J ^ 


^ T(a)da 


Then J can be reduced to the following Parseval relation 


(2tt) 


e ia r x(a)T (a)da 


where 


r = y - y 


(A. 5) 
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Now 


~ I -ia*x e ' ' 

T(a) = I e “T 


dx 


= 2tt I e 


TT 

/ 


a 


•iarcos(j)gi^^^^ J r 
o 


oo 

■s f 


- Xt 

e sinar dr 


Thus 


4tt 


4it 


— 5 7 ~ — 5 

+ X a + iio/v 


" ■ ” J 


. ->■ ->■ 

da 

“2T? 

to /V 


oo 

■t/ 


da 

a^ + co^/v^ 


IT 


27T r sin+d* 


03 

•/ 


g sin gr da 
g^ + to^/v"^ 


47^ 

(ojr/v) 


nr 
■ vir 


/ to ^ 


(A. 6) 


(A. 7) 
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APPENDIX B 


ACOUSTIC RADIATION FROM A VIBRATING SURFACE 
Consider the equations 


1 


o 



div v' 


0 


^ + grad h 


0 



f 


1 9^h' 



Let 


and 


Then 


and 


ll(x) = P_(h'(x.O") - h'(x,0'*')) 


h' ^ 


k = w/a. 


V^h' + k\' = 0 


Also 


h' = 


3z 


I 


Q(y) 


-ik|x-y I 
x-y I 


dy 


Lim h' = Q(x) = 
z->0± 



T dy 

R-^ i 


r 3 . (^2 ^ ^ 2 ) 3/2 


dy = pdpdG 


Lim h' = 
z->0± 


- 2 IT Q(x) 


£ 

/ 


pz dp 


( ' 1 , ■2V172 

(z +p ) 


= + 2tt Q(x) sgn z 


(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.IO) 


(B.ll) 
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and 


or 


Finally 


h'~ - h’+ = 4 tt q(x) = Mi) 

Poo 


Q(X) 


Ux) 



For |x| -> oo we get 


(B.12) 


(B.13) 


(B.14) 


where 


h' = 


1 , _9_ e 

47TP ’ 9z 




0 ik5c»y -> 

^(y) e dy 


X = 


‘x| 


r = Ixl 


Or 


h' 


ik 


z „-ikr 


/ 




(B.15) 


(B.16) 


(B.17) 


(B.18) 
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Also in the far field on a sphere of radius r 


3h' 

3r 




and 


i ^ 

(jj 9r 


3h' 

3r 


ik h' 


V. 


n 



h' 


a 


o 



(Re h 



; 2 2 

(Re h' = (h^cosojt - h^sinoot) 


= hl^ cos^wt+hl^ sin^oot 
^ J 


and so 


1 

I 


h' 


2 


h^ 



(B.19) 


(B.20) 


(B.21) 


(B.22) 


(B.23) 


(B.24) 


(B.25) 
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Ih' 

2 a. 


o 


■20 2 2 
J2tt p„a 


z2 


oo°'0 


5,(y)dy 


^''(y')dy' 


e^kx- (y-y ’) 
(B.26) 


sin0 = — 
r 


(B.27) 



Referring to the above sketch 


X - i cos6cos({) + j cosOsin^) + ^ sin6 


COS0 (Xg) + k sin0 


(B.28) 


X = X COS0 + j sincf) 


Thus 


.|2 


k^ sin^0 


■^a~ ~ 2 2 2" 

° 327T p^a^ r 


J" A(y)dy J", 


Uy)dy I 

S S 

(B.29) 
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Now 


Ea = 


d6 r COS0- 


-tt/2 


) )l(y)dy ) Jl“(y')ciy’ • G 

"^Piar, J J 


(B.30) 


32tt p^ao 


2tt tt/2 

- J do J d0 sln20coseei'>^‘=°=«^^s-<?-y'> (B.31) 


-tt/ 2 


For compact surfaces 


G = 2n I sin ecose 
-■ul2 
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